Charged Black Holes in Phantom Cosmology 
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Abstract 



In the classical relativistic regime, the accretion of phantom-like dark energy onto a 
■ stationary black hole reduces the mass of black hole. Here we have investigated the accretion 

CNI , of phantom energy onto a stationary charged black hole and have determined the condition 

under which this accretion is possible. This condition restricts the mass to charge ratio in a 
narrow limit. This condition also challenges the validity of the cosmic censorship conjecture 
, since a naked singularity is eventually produced as magnitude of charge increases compared 

' to mass of black hole. 

OO ■ 

O 

Keywords: Accretion; Black Hole; Phantom Energy. 

>< 

H , 

. P. , 1 Introduction 

Accelerated expansion of the universe has been observed and confirmed by myriad of sources 
including analysis of cosmic microwave background radiation [1], large scale structure [2] and 
supernovae SNe la data [3l H]. This expansion is supposedly driven by exotic vacuum energy 
having p > and p < or equivalently p = —p (or ui = —1), dominating the observable universe. 
Observations of WMAP data suggest that its magnitude is more than 70% of the total energy 
density of the universe [5]. Among other forms of exotic energies (e.g quintessence, cosmological 
constant, k-essence, hessence etc), the phantom energy with uo < —1 exhibits similar behavior 
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on large cosmic scale. The genesis of phantom energy is not clear but it violates the null and 
weak energy conditions. As these conditions are the weaker one, they ensure that the stronger 
conditions (i.e. strong and dominant) will be violated automatically [H [71 |8]. These energy 
conditions guarantee the positive definiteness of the energy densities and pressure densities of 
all the matter content in the universe. Recent observational data constrain the range of dark 
energy by —1.38 < uj < —0.82 at 95% confidence level [S]. Rather the supernovae data favor an 
evolving uj{z) varying from quintessence (w > —1) to phantom regime(a; < —1) [10]. Further, 
the extrapolation of WMAP data is best fitted with the notion of phantom energy [11] 

The energy density and the pressure of the phantom energy can be represented by the mini- 
mally coupled spatially homogeneous and time dependent scalar field (p having negative kinetic 
energy term given by 

p = -Y + n'A), p = -^-vi<i>). (1) 

Here V{(j)) is the scalar potential and dot over (p represents the derivative with respect to time 
parameter t. Note that if the kinetic term in Eq. (1) is positive then it gives usual dark energy 
with satisfies all the energy conditions. The above parameters p and p are related to the Hubble 
parameter H as 

H' = ^{-<P' + 2V), (2) 

where H(t) = a/a and a{t) is the scale factor which arises in the Friedmann-Robertson- Walker 
spacetime. From Eq. (2), we require the potential V{(j)) to be positive. It is argued by using 
scalar field models of phantom energy, that it can behave as a long range repulsive force [12j . 
The phantom energy possesses some peculiar properties unlike normal matter e.g. (1) its energy 
density p{t) increases with the expansion of the universe, (2) it ensures the existence and stability 
of traversable worm holes in the universe [13l [HI [151 [13 [H] > (3) also self-gravitating, static and 
spherically symmetric phantom scalar fields with arbitrary potentials can generate a stable 
configuration of a regular black hole or apparently non-singular black hole which inherently 
possesses exactly Schwarzschild-like causal structure but the singularity is replaced by a de Sitter 
infinity, thereby generating an asymptotically de Sitter expansion beyond the black hole horizon 
|181 119j. (4) due to strong negative pressure the phantom energy can disrupt all gravitationally 
bound structures i.e from galactic clusters to the gravitationally collapsed objects including 
black holes [201 [211 [22l [231 HH [25] , (5) it can produce infinite expansion of the universe in a finite 
time thus causing the 'big rip' (i.e. a state when a{t),p{t) ^ oo for t < cxd) [6, llj. 

The big rip is characterized by a future singularity implying a finite age of the universe. 
It has been proposed that this future singularity can be avoided if the phantom energy is 
interacting with the dark matter [261 [27] . The interaction of phantom energy and dark matter 
leads to stable attractor solutions at late times and the big rip is avoided in the parameter 
space |28) . Also, it was argued that there are certain classes of unified dark energy models 
stable against perturbations, in which cosmic dooms day can be avoided [29]. Moreover in 
scalar tensor theories, quantum gravity effects may prevent (or, at least, delay or soften) the 
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cosmic doomsday catastrophe associated with the phantom \30\ \31\ . Also in Gauss Bonnet 
gravity theory and loop quantum cosmology, the big rip occurrence is also avoided [SH [33] . 
In order to avoid the big rip with phantom matter, it is sufficient to have a phantom scalar 
field with a potential bounded above by some positive constant [Ml- It is also suggested that 
phantom dark energy with uj < —1 can effectively ameliorate the coincidence problem (i.e. why 
does the observable universe begin the accelerated expansion so recently and that why are we 
living in an epoch in which the dark energy and the matter energy density are comparable?) 
\35\ I36j. In another model using vector like dark energy with a background of perfect fluid, it is 
demonstrated that the cosmic coincidence problem is fairly solved |37| . 

The fate of a stationary uncharged black hole in the phantom energy dominated universe was 
investigated by Babichev et al [20] . The phantom energy was assumed to be a perfect fluid. 
The phantom energy was allowed to fall onto the black hole horizon only in the radial direction. 
It was concluded that black hole will lose mass steadily due to phantom energy accretion and 
disappear near the big rip. We here adopt their procedure for a static, stationary and charged 
black hole. The gravitational units are chosen for this work. 

The paper is organized as follows: In the second section, we have explained the relativistic 
model of accretion onto a charged black hole and obtained the black hole mass loss rate. In 
third section, we have determined the critical points of accretion model and have analyzed the 
dynamics about these points. Finally we conclude our paper. 



2 Accretion onto Charged Black Hole 

We consider a static and spherically symmetric black hole of mass M having electric charge e, 
so-called Reissner-Nordstrom (RN), specified by the line element 

ds'^ = f{r)dt^ - f{r)-^dr^ - r^{de^ + sin^ 9df^), (3) 

where 

If > then the metric is non-singular everywhere except at the curvature or the irremovable 
singularity at r = 0. Also if < then the function /(r) has two real roots given by 

rh± = M± - e2. (5) 

These roots physically represent the apparent horizons of the RN black hole. The two horizons 
are termed the inner rh- and the outer r/j+. The outer horizon is effectively called the event 
horizon while the inner one is called the cauchy horizon of the black hole. The metric (3) is 
then regular in the regions specified by the inequalities: oo > r > rh^, rh+ > r > r^- and 
r/i_ > r > 0. Note that if = M^, then it represents an extreme RN black hole while if 
> M^, it yields a naked singularity at r = [38\ I39j . 
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The phantom energy is assumed to be a perfect fluid specified by the stress energy tensor 



(6) 



Here p is the pressure and p is the energy density of the phantom energy. Also = {u''{r), u''(r), 0, 0) 
is the four velocity of the phantom fluid which satisfies the normalization condition u^n^ = — 1. 
We assume that the in-falling phantom fluid does not disturb the global spherical symmetry of 
the black hole. Further the energy-momentum conservation T.'^'^ = = T!'^ gives 



2M 



ur^M-^{p + p)\ 1 + _ + = (7i 



(7) 



where = u = dr / ds \s the radial component of the velocity four vector and Ci is a constant 
of integration. For inward flow, we will take n < 0. Moreover, the second constant of motion 
is obtained by projecting the energy conservation equation onto the velocity four vector as 
u^T^'J' = 0, which yields 



ur'^M ^ exp 



ph 



dp' 



P' + pip') 



-A. 



(8) 



Here ^ is a constant of integration. Above ph and poo are the energy densities of the phantom 
energy at the horizon and at infinity respectively. From Eqs. (7) and (8) we have 



, /, 2M e2 
[p + p)\ll — + ^ + u exp 



Ph 



dp' 



p' + p{p') 



C2, 



(9) 



where C2 = —Ci/A = p^o + p{poo)- In order to calculate the rate of change of mass of black 
hole we integrate the flux of the fluid over the entire cross-section of the event horizon as 



M 



T[dS, 



(10) 



where T[ determines the momentum density in the radial direction and dS = ^/—gd9dip is the 
surface element of the horizon, where g is the determinant of the metric. From Eqs. (7 - 10), 
we get 

A7rAM\p^{t)+Poo{t)), (11) 



dM 
IT 



which clearly demonstrates that mass of black hole decreases if p^o + Poo < 0. Note that Eq. 
(11) can be solved for any equation of state of the form p = p(p) or in particular p = ujp. In 
general, Eq. (11) holds for all p and p violating the dominant energy condition, thus we can 
write Ha iSI 



^=A7,AM\p{t)+p{t)). 



(12) 



In the astrophysical context, the mass of black hole is a dynamic quantity. The mass increases 
by the accretion of matter and can decrease by the accretion of the phantom energy. Since we 
are not incorporating matter in our model, the mass of black hole will decrease correspondingly. 
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3 Critical Accretion 



We are interested only in those solutions that pass through the critical point as these correspond 
to the material falling into the black hole with monotonically increasing speed. The falling fluid 
can exhibit variety of behaviors near the critical point of accretion, close to the compact object. 
For instance, for a given critical point r = rc, we have the following possibilities [IQ]: (a) = 
at r = Tc, — > as r ^ cxD, -u^ < for r > r^, and v? > for r < Vc- Thus for large distance, 
the speed of flow becomes negligible (subsonic), at the critical point it is sonic, while the flow 
becomes supersonic for very small r. Other solutions for the flow near are not of much interest 
due to their impracticality, like (b) u'^ < for all values of r and (c) > for all values of r. 
Solutions (b) and (c) are not realistic since they describe both subsonic and super-sonic flows 



for all r. Similarly, (d) 



for all values oi r > r^. and (e) = for all values of r < rc. 



Last two solutions are also useless since they give same value of speed at a given r. Hence from 
this discussion, we see that solution (a) is the only physically motivated, near the critical point. 
To determine the critical points of accretion we shall adopt the procedure as specified in Michel 
The equation of mass flux J.^ = gives 



pur = ki, 



(13) 



where ki is constant of integration. Dividing and then squaring Eqs. (7) and (13) give 

2 / ofl/f „2 \ \ 2 



P + P 
P 



2M n 
1 + ^ + 



k 



(14) 



Here C3 is a positive constant. Differentiation of Eqs. (13) and (14) and then elimination of dp 
gives 



du 
u 



2V^ 



M 
r 
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1 



or 



2M 



du 

dr 



+ h + u^ 



+ 



dr 



1 



2M 



+ h + u^ 



0, 



2^2 



M _ e2 

r 



N 

15' 



where 



y2 ^ d\n{p + p) _ ^ 
dlap 



(15) 



(16) 



(17) 



We have assumed that the flow is smooth at all points of spacetime, however if at any point the 
denominator D vanishes then the numerator N must also vanish at that point. Mathematically 
this point is called the critical point of the flow [H]. Equating the denominator D and numerator 
to zero, we can get the so-called critical point conditions given by 



Mrr 



2r2 



(18) 
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= .J , (19) 



and 

^2 _ Mr, - 

2r2 - 3Mrc + 

Note that by choosing e = in the above equations, we can retrieve the results for the accretion 
of fluid onto a Schwarzschild black hole [H]. All the quantities with subscript c are defined at 
the critical point correspondingly. Physically, the critical points represent the sonic point of the 
flow i.e. the point where the speed of flow becomes equal to the speed of sound, = or 
the corresponding Mach number = 1. This transition may occur from the initial subsonic 
to the supersonic or trans-sonic speeds. For any spherically symmetric spacetime, a surface 
where every point is a sonic point is called a sound horizon which itself will be spherical. Any 
perturbation or disturbance generated in the flow inside the sound horizon (r < r,) is eventually 
pulled towards the black hole singularity and hence cannot escape to infinity. 

It can be seen that the speed of sound (squared) = dp/ dp has no physical meaning if 
the EoS parameter w < (in p = ojp). Thus it will apparently make the exotic cosmic fluids 
like the cosmological constant, quintessence and the phantom energy unstable that can not be 
accreted onto the black hole. In order to avoid this problem, Babichev et al [l5] introduced a 
non-homogeneous linear equation of state (nEoS) given by p = a{p — Po), where the constants a 
and po are free parameters. The nEoS can describe both hydrodynamically stable {a > 0) and 
unstable {a < 0) fluids. The parameter uj is related to the nEoS as w = a{p — po) / p- Notice that 
a; < corresponds to a > and p > po, thus making the phantom energy as hydrodynamically 
stable fluid. Therefore the speed of sound Cs is now well-deflned with the nEoS for the phantom 
energy. Hence, the phantom energy can fall onto the RN black hole and can reduce the black 
hole mass. Since phantom energy reduces only mass and not charge, a stage is reached when 
the the cosmic censorship conjecture becomes violated i.e. e > m, the so-called emergence of a 
naked singularity. 

Now physically acceptable solution of Eq. (16) is obtained if > and > 0, hence we get 

- 3Mrc > 0, (20) 

and 

Mrc - > 0. (21) 

Eq. (20) can be factorized as 

2r2 - 3Mrc + = (r^ - r,+ )(re - r,_) > 0, (22) 

where 

rc± = ^(3M ± V9M2 _ 8e2), (23) 

which are positive satisfying rc+ > Vc- > 0. In general, for e < m, the inner critical point will 
lie between r/j_ < rc_ < r/j+, while the outer one will satisfy rc+ > rfi+- It is obvious that these 
roots will be real valued if 9M2 — 8e'^ > or 

^4- 
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These roots physically represent the locations of the critical or sonic points of the flow near the 
black hole. Notice that both mass and charge have same dimension of length, therefore all the 
inequalities here and below represent dimensionless ratios. Prom Eq. (22), we can see that these 
critical points specify two regions for the flow: (1) Tc > rc+ or (2) < rc < Tc-. We shall now 
solve Eq. (19) using (21) and then deduce a condition for the black hole mass and charge. 

To get solutions about the critical points, we substitute rc± in Eq. (21). For rc+, Eq. (21) 
gives 

MV9M2 - 8e2 > 4e^ - 3M^ (25) 

which is satisfled if 

^ < 1, (26) 

and 

4 

A comparison of inequalities (24), (26) and (27) imply 

8 M2 4 

Thus accretion through rc+ is possible if the above inequality (28) is satisfied. It encompasses the 
two types of black holes in itself: regular and the extreme RN black hole. Interestingly, the naked 
singularity also falls within the prescribed limits. Thus for all these spacetimes, the accretion 
is allowed through the critical point rc+. We stress here that using e = in the inequality (28) 
to retrieve same condition for the Schwarzschild black hole can be misleading. The inequality is 
deduced using the outer apparent horizon and a critical point. Since Schwarzschild black hole 
(e — > 0) possesses unique horizon and the critical point, the above inequality cannot be reduced 
for an uncharged black hole. 

Now we consider case (2) when < rc < rc_. Substitution of rc- in Eq. (21) gives 



MV'9M2 - 8e2 < 3M^ - 4e^. (29) 

If - 46^ < then Eq. (29) does not yield any solution. So we need SM^ - 46^ > which 
yields 

^ > 3- PO) 



Further inequality (29) is satisfied if 



M2 

< 1. (31) 



e2 



Since Eqs. (30) and (31) are mutually inconsistent, there is no solution for rc in case (2). Thus 
accretion is not possible through rc- . 

Since the mass of black hole is decreasing by the accretion of phantom energy (see Eq. 12), 
it implies that at least one critical point must exist for the fiuid fiow, which is specified by rc+. 
This critical point yields the mass to charge ratio of the black hole in the range specified by (28) 
which allows that accretion onto all charged spherically symmetric black holes. 
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4 Conclusion 



We have analyzed the effects of accretion of phantom energy onto a charged black hole. The 
analysis is performed using two critical points rc±. It turns out that accretion is possible 
only through rc+ which yields a constraint on the mass to charge ratio given by Eq. (28). 
This expression incorporates both extremal and non-extremal black holes. Thus all charged 
black holes will diminish near the big rip. Apparently this condition predicts the existence of 
large charges onto black holes, although astrophysically no such evidence has been successfully 
deduced from the observations. In theory, the existence of large charges onto black holes is 
consistently deduced by the general theory of relativity. It needs to be stressed that there is 
no analogous condition for the Schwarzschild black hole (e = 0). This analysis can be extended 
for a rotating charged black hole (so-called Kerr-Neumann black hole) to get a deeper insight of 
the accretion process. This work also serves as the generalization of Michel [H] in terms of the 
accretion of phantom dark energy onto a charged black hole. 
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